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Abstract 

The history of the transition from a neutral intergalactic medium (IGM) to one that 
■ is almost fully ionized can reveal the character of cosmological ionizing sources and set 

important constraints on the stellar birthrate at high redshifts. The hydrogen component 
in a highly inhomogeneous universe is completely reionized when the number of photons 
emitted above 1 ryd in one recombination time equals the mean number of hydrogen 
atoms. If stellar sources are responsible for photoionizing the IGM at z = 5, the rate of 
star formation at this epoch must be comparable or greater than the one inferred from 
optical observations of galaxies at z ~ 3, and the mean metallicity per baryon in the 
universe <; 1/500 solar. In hierarchical clustering scenarios, high-z dwarfs (i.e. an early 
generation of stars in dark matter halos with circular velocities v c i rc ~ 50 kms" 1 ) are 
expected to be one of the main source of UV photons and heavy elements at early epochs. 
They would be very numerous, > 0.2 arcsec -2 , and faint, Tab <^ 29.5 mag: their detection 
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may have to wait for the Next Generation Space Telescope. 



1 Introduction 



What keeps the universe ionized at z = 5? The existence of a filamentary, low-density inter- 
galactic medium (IGM), which contains the bulk of the hydrogen and helium in the universe, 
is predicted as a product of primordial nucleosynthesis || and of hierarchical models of grav- 
itational instability with "cold dark matter" (CDM) jp, |47_|], [17]. The application of the 



Gunn-Peterson |10| constraint on the amount of smoothly distributed neutral material along 
the line of sight to distant objects requires the hydrogen component of the diffuse IGM to have 
been highly ionized by z ~ 5 f37f , and the helium component by z ~ 2.5 ||. The plethora of 



discrete absorption systems which give origin to the Lya forest in the spectra of background 
quasars are also inferred to be strongly photoionized. From QSO absorption studies we also 
know that neutral hydrogen accounts for only a small fraction, ~ 10%, of the nucleosynthetic 
baryons at early epochs fl9fl . It thus appears that substantial sources of ultraviolet photons 
were present at z > 5, perhaps low- luminosity quasars (TB| or a first generation of stars in 



virialized dark matter halos with T vir ~ 10 4 — 10 5 K ||, [32], |[T^], p9 |: early star formation 
provides a possible explanation for the widespread existence of heavy elements in the IGM J7| . 
More in general, establishing the character of cosmological ionizing sources is an efficient way to 
constrain competing models for structure formation in the universe, and to study the collapse 
and cooling of small mass objects at early epochs. 

In this talk I will focus on the candidate sources of photoionization at early times and on 
the time-dependent reionization problem, i.e. on the history of the transition from a neutral 
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Figure 1: Comoving space density of bright QSOs as a function of redshift. The data points 



with error bars are taken from [|T^] (filled dots), |45[ (filled squares), |36| (crosses), and [fig 
(filled pentagon). The points have been normalized to the z = 2.5 space density of quasars 
with Mb < —26 (Mb < —27 in the case of |18|). The empty triangles show the space density 
(normalized to the peak) of the Parkes flat-spectrum radio-loud quasars with P > 7.2 x 10 26 
W Hz" 1 sr" 



IGM to one that is almost fully ionized. Throughout this paper I will adopt an Einstein-de 
Sitter universe (go = 0.5) with H = 50h 50 kms -1 Mpc -1 . The ideas described below have been 



developed in collaboration with F. Haardt and M. J. Rees |25 . 



2 Sources of Ionizing Radiation at High Redshifts 

2.1 Quasars 

The existence of a decline in the space density of bright quasars at redshifts beyond ~ 3 was first 
suggested by and has been since then the subject of a long-standing debate. In recent years, 
several optical surveys have consistently provided new evidence for a turnover in the QSO counts 
|T~5j| , f45| , f36j, |18[]. The interpretation of the drop-off observed in optically selected samples 
is equivocal, however, because of the possible bias introduced by dust obscuration arising from 



intervening systems | 33fl . Radio emission, on the other hand, is unaffected by dust, and it has 



recently been shown |39| that the space density of radio-loud quasars also decreases strongly 
for z > 3, demonstrating that the turnover is indeed real and that dust along the line of sight 
has a minimal effect on optically-selected QSOs (Figure 1). 

The QSO emission rate of hydrogen ionizing photons per unit comoving volume can be 



written as 
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where a s is the slope of the quasar spectral energy distribution, and the best-fit parameters 
are = 1.9, ( = 2.58, and £ = 3.16. It is important to notice that the procedure adopted 
to derive this quantity implies a large correction for incompleteness at high-z. With a fit to 
the quasar luminosity function (LF) which goes as <f)(L) oc L -L64 at the faint end [j3~4] , the 



contribution to the emissivity converges rather slowly, as L 36 . At z — 4, for example, the 
blue magnitude at the break of the LF is M* « —25.4, comparable or slightly fainter than the 
limits of current high- ,2 QSO surveys. A large fraction, about 90% at z = 4 and even higher 
at earlier epochs, of the ionizing emissivity in our model is therefore produced by quasars that 
have not been actually observed, and are assumed to be present based on an extrapolation 
from lower redshifts. The value of Mq obtained by including the contribution from observed 
quasars only would be much smaller at high redshifts than shown in Figure 2. While it is also 
possible that an excess of low-luminosity QSOs, relative to the best-fit LF, could actually boost 
the estimated ionizing emissivity at early epochs, the observed lack of red, unresolved faint 
objects in the Hubble Deep Field seems to argue against models where the quasar LF steepens 



significantly with lookback time |14 . 
2.2 Star-forming Galaxies 

Galaxies with ongoing star-formation are another obvious source of Lyman continuum photons. 
The tremendous progress in our understanding of faint galaxy data made possible by the recent 
identification of star-forming galaxies at 2 < z < 4 in ground-based surveys K2j and in the 



Hubble Deep Field (HDF) jy]], [p4 |, |22| has provided new clues to the long-standing issue 
of whether galaxies at high redshifts can provide a significant contribution to the ionizing 
background flux. Since the rest-frame UV continuum at 1500 A (redshifted into the visible band 
for a source at z ~ 3) is dominated by the same short-lived, massive stars which are responsible 
for the emission of photons shortward of the Lyman edge, the needed conversion factor, about 
one ionizing photon every 10 photons at 1500 A, is fairly insensitive to the assumed IMF and 
is independent of the galaxy history for t ^> 10 7 3 yr. A composite ultraviolet luminosity 
function of Lyman-break galaxies at z ~ 3 has been recently derived by ||. It is based on 
spectroscopically (about 375 objects) and photometrically selected galaxies from the ground- 
based and HDF samples, and spans about a factor 50 in luminosity from the faint to the bright 
end. Because of the uncertanties that still remain in the rescaling of the HDF data points 
to the ground-based data, the comoving luminosity density at 1500 A is estimated to vary 
within the range 1.6 to 3.5 x 10 26 ergss -1 Hz -1 Mpc -3 . The "best guess" Schechter fit gives 
log (0*/Gpc -3 ) = 6.1, a = 1.38, and MX B (1500) = -20.95 @, the magnitude at the "break" 
corresponding to a star-formation rate slightly in excess of 10M Q yr -1 (Salpeter IMF). 

Figure 3 shows the Lyman continuum luminosity function of galaxies at z ~ 3 (at all ages 
^0.1 Gyr one has L(1500) / L(912) « 6 for a Salpeter mass function and constant star formation 
rate 0), compared to the distribution of QSO luminosities at the same redshift. The comoving 
ionizing emissivity due to Lyman-break galaxies is 4.2 ±1.5 x 10 25 ergss -1 Hz -1 Mpc -3 , between 
2 and 4 times higher than the estimated quasar contribution at z — 3. 

This number neglects any correction for dust extinction and intrinsic H i absorption. While 



it has been pointed out by many authors ||28[| , [ 55] , || that the colors of Lyman-break galaxies 
are redder than expected in the case of dust-free star-forming objects, the prescription for a 
"correct" de-reddening is still unknown at present (note that redder spectra may also results 




redshift 



Figure 2: Comoving emission rate of hydrogen Lyman continuum photons (solid line) from 
QSOs, compared with the minimum rate (dashed line) which is needed to fully ionize a fast 
recombining (with clumping factor C = 30) Einstein-de Sitter universe with h 50 = 1 and 
Qb = 0.08. Models based on photoionization by quasar sources appear to fall short at z = 5. See 
for details on the assumed quasar luminosity function and spectral energy distribution. The 
data points with error bars show the estimated contribution of star-forming galaxies at z ~ 3 
and, with significantly larger uncertainties, at z ~ 4. The fraction of Lyman continuum photons 
which escapes the galaxy H i layers into the intergalactic medium is taken to be / e sc = 0.5. 



from an aging population or an IMF which is rather poor in massive stars). A Salpeter IMF, 
Eb-y = 0.1 model with SMC- type dust in a foreground screen, for example, has been found 



to reproduce quite well the rest-frame ultraviolet colors of the HDF "UV dropouts" ||26|| . In 
this model the color excess S912-1500 = 1.64Eb-v is rather small and can be safely neglected 
in correcting from observed rest-frame far-UV to the Lyman edge. For typical dust-to-gas 
ratios, however, it is the H 1 associated with dust that would greatly reduce the flux of Lyman 
continuum photons able to escape into the intergalactic space. The data points plotted in 
Figure 2 assume a value of / e sc = 0.5 for the unknown fraction of ionizing photons which 



escapes the galaxy H 1 layers into the intergalactic medium [27|. The possible existence of 
a numerous population of galaxies below the detection threshold, i.e. having star formation 
rates < 0.5M Q yr~ 1 , with a space density well in excess of that predicted by extrapolating to 
faint magnitudes the a = 1.38 best-fit Schechter function, will be discussed below. The LF 
of Lyman-break galaxies at z ~?Z 4 is highly uncertain. An analysis of the 5-band dropouts in 
the HDF - candidate star-forming objects at 3.5 < z < 4.5 - seems to imply a decrease in 
the comoving UV galaxy emissivity by about a factor of 2.5 in the interval 2.75 < z < 4 [23], 
|[26|| . In this sense star-forming galaxies with SFR in excess of 0.5 M yr" 1 may have a negative 
evolution with lookback time similar to the one observed in bright QSOs, but the error bars are 
still rather large. Adopting a L(1500) to L(912) conversion factor of 6, we estimate a comoving 
ionizing emissivity of 1.7 ± 1.1 x 10 25 / e scergss _1 Hz -1 Mpc -3 at z ~ 4. One should note that, 




Figure 3: The 912 A luminosity function of galaxies at z fa 3 (solid line), compared to the dis- 
tribution of QSO luminosities at the same redshift (dashed line). The latter has been derived 
assuming a spectral slope of ot s = 0.5. The former assumes a Salpeter IMF with constant con- 
stant star formation rate (age=l Gyr): Mab(912 A)= —19 corresponds to a rate of 13M Q yr _1 . 
The solid and dashed lines represent functional fits to the data points, and the dotted lines 
their extrapolation. 

while a population of highly reddened galaxies at high redshifts would be missed by the dropout 
color technique (which isolates sources that have blue colors in the optical and a sharp drop in 
the rest-frame UV), it seems unlikely that very dusty objects (with / esc <C 1) would contribute 
in any significant manner to the ionizing metagalactic flux. 



3 Reionization of the Universe 

In inhomogeneous reionization scenarios, the history of the transition from a neutral IGM to 
one that is almost fully ionized can be statistically described by the evolution with redshift 
of the volume filling factor or porosity Q(z) of H n , He n , and He in regions. The radiation 
emitted by spatially clustered stellar-like and quasar-like sources - the number densities and 
luminosities of which may change rapidly as a function of redshift - coupled with absorption 
processes in a medium with a time- varying clumping factor, all determine the complex topology 
of neutral and ionized zones in the universe. When Q << 1 and the radiation sources are 
randomly distributed, the ionized regions are spatially isolated, every UV photon is absorbed 
somewhere in the IGM, and the ionization process cannot be described as due to a statistically 
homogeneous radiation field . As Q grows, the crossing of ionization fronts becomes more and 
more common, and the neutral phase shrinks in size until the reionization process is completed 
at the "overlap" epoch, when every point in space is exposed to Lyman continuum radiation 
and Q — 1. 



When an isolated point source of ionizing radiation turns on, the ionized volume initially 
grows in size at a rate fixed by the emission of UV photons, and an ionization front separating 
the H ii and H i regions propagates into the neutral gas. Most photons travel freely in the 
ionized bubble, and are absorbed in a transition layer. The evolution of an expanding H n region 
is governed by the equation 

( 2 ) 

at Tift t rec 

381, where Vj is the proper volume of the ionized zone, N ion is the number of ionizing photons 



emitted by the central source per unit time, H is the Hubble constant, is the mean hydrogen 
density of the expanding IGM, n H (0) = 1.7 x 10~ 7 (n^/0.08), and 



i iec = [(l + 2 X )n H a B Cr = 0.3Gyr(^ VllfV'ctf (3) 



is the volume-averaged gas recombination timescale. Here olb is the recombination coefficient 
to the excited states of hydrogen, x t ne helium to hydrogen cosmic abundance ratio, C = 
( n Hii)/^Hii i s the ionized hydrogen clumping factor, [] and a gas temperature of 10 4 K has been 
assumed. Clumps which are dense and thick enough to be self-shielded from UV radiation will 
stay neutral and will not contribute to the recombination rate. An empirical determination 
of the dumpiness of the IGM at high redshifts is hampered by our poor knowledge of the 
ionizing background intensity and the typical size and geometry of the absorbers. Numerical 
N-body/hydrodynamics simulations of structure formation in the IGM within the framework 
of CDM dominated cosmologies have recently provided a definite picture for the origin of 
intervening absorption systems, one of an interconnected network of sheets and filaments, with 
virialized systems located at their points of intersection. In the simulations of |J, for example, 
the clumping factor rises above unity when the collapsed fraction of baryons becomes non 
negligible, i.e. z < 20, and grows to C > 10 (40) at z ~ 8 (5) (because of finite resolution effects, 
numerical simulations will actually underestimate clumping): the recombination timescale is 
much shorter than that for a uniform IGM, and always shorter than the expansion time. 

When t rec <C t, the growth of the H n region is slowed down by recombinations in the highly 
inhomogeneous IGM, and its evolution can be decoupled from the expansion of the universe. 
Just like in the static case, the ionized bubble will fill its time-varying Stromgren sphere after 
a few recombination timescales, 

Vi= AU^c (1 _ e _ t/froc) _ (4) 
n H 

One should point out that the use of a volume-averaged clumping factor in the recombination 
timescale is only justified when the size of the H n region is large compared to the scale of 
the clumping, so that the effect of many clumps (filaments) within the ionized volume can be 
averaged over. This will be a good approximation either at late epochs, when the IGM is less 
dense and H n zones have had time to grow, or at earlier epochs if the ionized bubbles are 
produced by very luminous sources like quasars or the stars within halos collapsing from very 
high-a peaks. The mean free path between absorbers having neutral columns > iVm is 0.8 
Mpc h^g 1 [(1 + z)/6)] -4 ' 5 (./Vhi/10 15 cm" 2 ) 0,5 : it is only on scales greater than this value that the 
clumping can be averaged over. On smaller scales underdense regions are ionized first, and only 
afterwards the UV photons start to gradually penetrate into the higher density gas. 

1 This may be somewhat lower than the total gas clumping factor if higher density regions are less ionized 



With these caveats in mind, equation approximately holds for every isolated source of 
ionizing photons in the IGM. The filling factor of H n regions in the universe, Qhii ; is then equal 
at any given instant t to the integral over cosmic time of the rate of ionizing photons emitted 
per hydrogen atom and unit cosmological volume by all radiation sources present at earlier 
epochs, Jq hi OQ (t')dt' /nn(t'), minus the rate of radiative recombinations, Jq Qmi(t')dt' /i Tec (t'). 
Differentiating one gets 

dQuu "ion Qhii 

dt 7Z H tree ' 

It is this simple differential equation - and its equivalent for expanding helium zones - that sta- 
tistically describes the transition from a neutral universe to a fully ionized one, independently, 
for a given emissivity, of the complex and possibly short-lived emission histories of individual 
radiation sources, e.g., on whether their comoving space density is constant or actually varies 
with cosmic time. In the case of a time-independent clumping factor, equation (|5|) has formal 
solution 

QhiiW = f dt' exp [-1- + £L) , (6) 



'H'H y ^rcc ^rec^ y 

with i rec oc t' 2 . At high redshifts, and for an IGM with C 1, one can expand around t to find 

Qffll(t) « —tree (7) 

The porosity of ionized bubbles is then approximately given by the number of ionizing photons 
emitted per hydrogen atom in one recombination time. In other words, because of hydrogen 
recombinations, only a fraction t iec /t (~ a few per cent at z — 5) of the photons emitted above 
1 ryd is actually used to ionize new IGM material. The universe is completely reionized when 
Q — 1, i.e. when 

^iorA-ec = ^H- (8) 

While this last expression has been derived assuming a constant comoving ionizing emissivity 
and a time-independent clumping factor, it is also valid in the case n ion and C do not vary 
rapidly over a timescale i Tec . 



4 Dwarfs at High Redshift? 

As t rec <C t at high redshifts, it is possible to compute at any given epoch a critical value for 
the photon emission rate per unit cosmological comoving volume, A/jon; independently of the 
(unknown) previous emission history of the universe: only rates above this value will provide 
enough UV photons to ionize the IGM by that epoch. One can then compare our determinations 
of A/ion to the estimated contribution from QSOs and star-forming galaxies. Equation @ can 
then be rewritten as 

The uncertainty on this critical rate is difficult to estimate, as it depends on the clumping 
factor of the IGM (scaled in the expression above to the value inferred at z = 5 from numerical 
simulations 0) and the nucleosynthesis constrained baryon density. A quick exploration of the 
available parameter space indicates that the uncertainty on A/ion could easily be of order ±0.2 in 
the log. The evolution of the critical rate as a function of redshift is plotted in Figure 2. While 
A/ion is comparable to the quasar contribution ai z > 3, there is some indication of a significant 



deficit of Lyman continuum photons at z — 5. For bright, massive galaxies to produce enough 
UV radiation at z — 5, their space density would have to be comparable to the one observed at 
z ~ 3, with most ionizing photons being able to escape freely from the regions of star formation 
into the IGM. This scenario may be in conflict with direct observations of local starbursts 
below the Lyman limit showing that at most a few percent of the stellar ionizing radiation 
produced by these luminous sources actually escapes into the IGM [|(| .f\ If, on the other hand, 
faint QSOs with (say) Mab = — 19 at rest-frame ultraviolet frequencies were to provide all the 
required ionizing flux, their comoving space density would be such (0.0015 Mpc -3 ) that about 
50 of them would expected in the HDF down to Jab = 27.2. At z > 5, they would appear very 
red in V — I as the Lya forest is shifted into the visible. This simple model can be ruled out, 
however, as there is only a handful (7) of sources in the HDF with (V — I)ab > 1-5 mag down 
to this magnitude limit. 

It is interesting to convert the derived value of A/i on into a "minimum" star formation rate 
per unit (comoving) volume, (hereafter we assume ^ b /2j = 0.08 and C = 30): 

P*{z) =tiUz) x lO- 53 " 1 ^ 1 « 0.013/^ (^) 3 Meyr^Mpc" 3 . (10) 

The conversion factor assumes a Salpeter IMF with solar metallicity [0. It can be understood 
by noting that, for each 1 M of stars formed, 8% goes into massive stars with M > 2OM 
that dominate the Lyman continuum luminosity of a stellar population. At the end of the 
C-burning phase, roughly half of the initial mass is converted into helium and carbon, with a 
mass fraction released as radiation of 0.007. About 25% of the energy radiated away goes into 
ionizing photons of mean energy 20 eV. For each 1 M & of stars formed every year, we then 
expect 

0.08 x 0.5 x 0.007 x 0.25 x M c 2 1 53 



10 53 phots" 1 (11) 
20 ev lyr 

to be emitted shortward of 1 ryd. Note that the star formation density given in equation fllOD 
is comparable with the value directly "observed" (i.e., uncorrected for dust reddening) at z ~ 3 

B I- 

The same massive stars that dominate the Lyman continuum flux also manufacture and 
return most of the metals to the ISM. In the approximation of instantaneous recycling, the rate 
of ejection of newly sinthesized heavy elements which is required to keep the universe ionized 
at redshift z is, from equation fllP]), 



p z (z)=y(l-R)p*(z)> 3.5xl0- 4 ( JLj f£ M yr" 1 Mpc" 3 , (12) 



where y is the net, IMF-averaged "yield" of returned metals, Z Q = 0.02, and R ~ 0.3 is the 
mass fraction of a generation of stars that is returned to the interstellar medium. At z — 5, 
and over a timescale of At = 0.5 Gyr corresponding to a formation redshift Zf = 10, such a 
rate would generate a mean metallicity per baryon in the universe of 

< z >~ 3H%: * om2 ufe ^ (13) 



more 



comparable with the level of enrichment observed in the Lya forest at z ~ 3 pOf : 
than 2% of the present-day stars would need to have formed by z ~ 5. It has been recently 
suggested p9l that a large number of low-mass galactic halos, expected to form at early times 
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Note that, at z = 3, Lyman-break galaxies would radiate more ionizing photons than QSOs for / csc > 30%. 



in hierarchical clustering models, might be responsible for photoionizing the IGM at these 
epochs. According to spherically-symmetric simulations f44jl , photoionization heating by the 



UV background flux that builds up after the overlapping epoch completely suppresses the 
cooling and collapse of gas inside the shallow potential wells of halos with circular velocities 
< 35 kms" 1 . Halos with circular speed v c \ rc = 50 kms -1 , corresponding in top-hat spherical 
collapse to a virial temperature T vir = 0.5^m p v^ iTC /k »s 10 5 K and halo mass M = 0.1v^ iTC /GH « 
4 x 10 9 [(1 + z) /6]~ 3 ^ 2 h^Q M , appear instead largely immune to this external feedback (but see 
[p0| ). In these systems rapid cooling by atomic hydrogen can then take place and a significant 
fraction, /fi&, of their total mass may be converted into stars over a timescale At comparable 
to the Hubble time (with / close to 1 if the efficiency of forming stars is high). If high-z 
dwarfs with star formation rates /Q^M/At ~ 0. 3/0.5 Ai^5 M yr _1 were actually responsible 
for keeping the universe ionized at z ~ 5, their comoving space density would have to be 

0.013/^ Me yr-Mpc- /JU/\- 3 , (14) 

0.6/MQyr- 1 ^0.25/ 1 ' V ; 

two hundred times larger than the space density of present-day galaxies brighter than L*(4400), 
and about five hundred times larger than that of Lyman-break objects at z ~ 3 with M < 
M AB (1500), i.e. with star formation rates in excess of 10 M yr _1 . Only a rather steep lumi- 
nosity function, with Schechter slope a ~ 2, would be consistent with such a large space density 
of faint dwarfs and, at the same time, with the paucity of brighter B- and V^-band dropouts 
observed in the HDF. The number density on the sky would be 0.2 arcsec -2 , corresponding 
to more than three thousands sources in the HDF. With a typical apparent magnitude at z = 5 
of Jab ~ 29.5 mag (assuming / = 0.5), these might be too faint to be detected by the HST, 
but within the range of the proposed Next Generation Space Telescope . 



These estimates are in reasonable agreement (although towards the low side) with the 



calculations of |29| , who used the observed metallicity of the Lya forest as the starting point of 
their investigation. A higher density of sources - which would therefore have to originate from 
lower amplitude peaks - would be required if the typical efficiency of star formation and/or 
the escape fraction of ionizing photons were low, (/, / esc ) < 1. In this case the dwarfs could 
still be detectable if a small fraction of the gas turned into stars in very short bursts (there 
would then be an extra parameter associated with their duty cycle, in addition to / esc and /). A 
reduction of the star formation rate in halos with low circular velocities (necessary in hierarchical 
cosmogonies to prevent too many baryons from turning into stars as soon as the first levels of 
the hierarchy collapse f46f ) may result from the heating and possible expulsion of the gas due 
to repeated supernova (SN) explosions after only a very small number of stars have formed. 
Recent numerical simulations p|] show, however, that metal-enriched material from SN ejecta 
is accelerated to velocities larger than the escape speed from such systems far more easily than 
the ambient gas. If the same population of dwarf galaxies that may keep the universe ionized at 
z w 5 were also responsible for polluting the IGM with heavy elements, it is interesting to ask 
whether these atoms could diffuse uniformly enough to account for the observations of weak but 
measurable C iv absorption lines associated with the Lya forest clouds pl| , [jjQfl . Our fiducial 



estimate of 0.1 sources per comoving Mpc 3 implies that, in order to pollute the entire IGM, the 
metal-enriched material would have to be expelled to typical distances of order 1 Mpc, which 
would require an ejection speed of about 200 [(1 + z)/6} 1 ^ 2 kms -1 . In fact, the heavy elements 
may be restricted to filaments within only (say) 20 proper kpc (100-150 comoving kpc) of the 
halos, and in this case it would be enough to accelerate the outflowing metals to velocities 
comparable to the escape velocity, rather than to the higher speeds associated with SN ejecta. 
The required diffusion of heavy elements would be a more serious constraint if the relevant 



galaxies were rarer and more luminous, as would happen if they originated from 3-<r peaks and 
star formation was postulated to occur with higher efficiency than in more typical peaks. 
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